Sharp interface limit of a phase-field model of crystal grains 
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We analyze a two-dimensional phase field model designed to describe the dynamics of crystalline 
grains. The phenomenological free energy is a functional of two order parameters. The first one 
reflects the orientational order while the second reflects the predominant local orientation of the 
crystal. We consider the gradient flow of this free energy. Solutions can be interpreted as ensembles 
of grains (in which the phase of the order parameter is approximately constant in space) separated 
by grain boundaries. We study the dynamics of the boundaries as well as the rotation of the grains. 
In the limit of the infinitely sharp interface, the normal velocity of the boundary is proportional to 
both its curvature and its energy. We obtain explicit formulas for the interfacial energy and mobility 
and study their behavior in the limit of a small misorientation. We calculate the rate of rotation of 
a grain in the sharp interface limit and find that it depends sensitively on the choice of the model. 
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I. INTRODUCTION 

The characterization and evolution of microstructure 
forms a cornerstone of materials science. In particular the 
grain structure of a polycrystalline material determines 
many of its properties. Recent efforts at modeling the 
evolution of grain boundaries have used a variety of ap- 
proaches 0] , H , H , H . Herein we focus on the recently 
introduced phase field model of Kobayashi, Warren and 
Carter (KWC)|| which is based on earlier attempts by 
the same authors (M and 0). 

The model is motivated by symmetry principles and 
has a surprisingly rich set of physical characteristics. In 
particular, KWC showed numerically that solutions to 
this model can be interpreted as a collection of grains. 
The velocity of the interface between the grains was 
found to be approximately proportional to the local cur- 
vature of the interface, but the grains were also able to 
rotate towards lower energy misorientations. In support 
of the notion of grain rotation, three independent unpub- 
lished molecular dynamics studies by M. Upmanyu and 
D. Srolovitz, S. R. Phillpot and D. Wolf, and S. Srivilla- 
puthur and J. W. Cahn, suggest that grain rotation will 
occur under certain circumstances. In addition, there is 
a substantial history (and debate) concerning the mech- 
anisms [U , M and observation of grain rotation |1Q, 1 1 



The KWC phase field model is challenging mathe- 
matically because of a singular term in the free energy. 
Kobayashi and Giga studied similar singular models 
and showed that there is a way to handle the singularity 
consistently. In this paper we will apply their method to 
the KWC model and show that its solutions can indeed 
be interpreted as grains. Our main goal is to analytically 
obtain the properties of a grain boundary as well as the 
rotation rate of a grain. We accomplish this task by con- 
sidering a distinguished limit of the model parameters in 
which the width of the boundary vanishes while its mea- 
surable characteristics remain finite and non-zero. The 
methodology of this this so called sharp interface limit is 
well established fl§, @, fy|. 

The organization of this paper is as follows. In Sec. II 
we introduce the order parameters, phenomenological 



free energy, and the gradient flow equations. We next 
perform a formal asymptotic expansion of the model in 
Sec. III. The zcroth order probelm of this expansion is 
discussed in Sec. IV, where we obtain the profile and the 



energy of the static flat boundary. In the following sec- 
tion (Sec. |v|) the first order asymptotics are examined. 
We are able to determine the velocity of a curved grain 
boundary. We find that this velocity is proportional to 
both the curvature and the interfacial energy and obtain 
an expression for the mobility of the interface. As noted 
above, grains can rotate in this model. We calculate the 
rate of rotation of a grain in Sec. VI. Finally, in Sec. VII 



we apply all of the results described herein to the simple 
case of a circular grain embedded in a matrix. In our 
conclusion we discuss further ideas and the implications 
of this phase field m odel for the problem of coarsening in 



polycrystals in Sec. VIII 



II. MODEL 



We model the evolution of a collection of nearly per- 
fect crystalline grains in two dimensionsvia a phase field 
model. First, we discuss order parameters which cap- 
ture the microscopic physics of grain boundaries. It is 
then possible to construct a phenomenological free energy 
which favors a perfect uniform crystal and supports sta- 
ble grain boundaries. Evolution of an ensemble of grains 
is then modeled via gradient flow of this free energy. 

Following ^| we develop two order parameters which 
capture the physics of grain boundaries. To distinguish 
grains of different orientations we introduce a continu- 
ously varying local orientation 8. Since the energy of 
crystal does not depend on 9 itself, the phenomenological 
free energy will be a functional of the gradients of only. 
The second order parameter 77 is used to differentiate the 
nearly perfect crystal in the interior of the grains from 
the disordered material in the grain boundary. It varies 
from perfect order ij = 1 to complete disorder rj = 0. 
Both order parameters are not conserved. 
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We analyze the free energy (based on KWC) 



- / dA 



or , l9 



f(r))+ g(ri)8\V0\ 



(1) 



where a, e and s are positive model parameters. Some 
readers may find the overall prefactor of 1/e suggestive. 
Subsequently we will examine the limit e — > 0. This pref- 
actor ensures that the surface energy of a grain boundary 
tends to a non-zero constant. 

Term by term, the above free energy deserves some 
discussion, although the reader interested in the full mo- 
tivation of this model is referred to ||. The first term 
describes the penalty for gradients in the order parameter 
(grain boundaries cost energy). The free energy density 
f(rj) is chosen to be a single well with the minimum at 
rj = 1 and /(l) = reflecting the fact that disordered 
material has higher free energy. The third and fourth 
terms are an expansion in |V#|, where the couplings g(r)) 
and h(rf) must be positive definite. KWC argued that the 
expansion in |V0| must begin at first order to assure the 
existence of stable grain boundaries. This term yields a 
singularity in the dynamic equations. 

Indeed, Ref. || omitted the |V6*| 2 term in their analyt- 
ical investigation of a stationary flat interface. This term 
was added in 0, for practical reasons, in order to solve 
the model numerically. We shall see that this term makes 
no qualitative difference in the properties of static grain 



boundaries. However, as shown in Sec. V A, the mobility 
of a grain boundary vanishes as h 1 / 2 when in the h — ► 
limit. Thus, the extra term introduced in KWC seems to 
be required for grain boundaries to move. 

Assuming relaxational dynamics for a non-conserved 
set of order parameters, we find the gradient flow equa- 
tions: 



di] 
dt 



8T_ 

5rj 



(2a) 



a 2 V 2 77 -/„- fl „s|V0| -h v -\V9\ 2 , 



P( V ,V0)Tgr 1 : 



89 

~dt 



- 2 V 2 7 

5JF_ 
'56 



(2b) 



h e 2 V9 + g s 



V6> 

W\ 



where we used the subscript rj to denote differentiation. 
The mobility functions P and Q must be positive def- 
inite, continuous at V# = 0, but are otherwise unre- 
stricted. The system (^) must be supplemented by ini- 
tial and boundary conditions and a rule that specifies the 
handling of the indctcrminancy (and sigular divergence) 
of the term V0/| V6>|. This particular type of the singular- 
ity is generally handled in a theory of extended gradients. 
Ref. proves that there is indeed a unique way to pre- 
scribe the value of the right hand side of Eq. (2b) when 
V9 = 0. While we do not wish to attempt to explain all 



of the details of Ref. |l2j , it is useful to summarize the 
ultimate conclusion of the mathematical analysis. Thus, 
let us define a collection {Gi} of distinct connected re- 
gions where V0 = 0. We shall refer to Gi as the interior 
of grain i. The essence of the method of dealing with the 
singularity is that the orientation 9 must remain uniform 
in space in each Gi. Therefore, the right hand side of 
Eq. ( pb| ) is chosen in to be uniform in Gi. This condi- 
tion, along with the requirement that 77, 9 and V(9 be 
continuous at the boundaries of Gi, uniquely determines 
the rotation rate of each grain (i.e. 89 /dt in each Gi) and 
the motion of its boundaries. 




FIG. 1: Two-dimensional simulation of Eqs. (g) with simple 
choices for the couplings. The orientation 8 is constant in the 
shaded grain regions. The color is used to represent the value 
of 6. Dark regions correspond to small values of r). 



Ref. Q verifies numerically that typical solutions of 
the gradient flow system indeed consist of a collection 
of connected regions Gi in which V# = and rj w 1 (grain 
interiors) separated by narrow regions (grain boundaries) 
in which the orientation 9 changes smoothly between the 
neighboring grains. Thus, the grains rotate (9 changes 
in time uniformly in each grain), and grain boundaries 
migrate (see Fig. |l] adopted from Ref. p6[). 



III. FORMAL ASYMPTOTICS 

Having elucidated the general properties of the model, 
we now examine the sharp- interface limit, in order to ex- 
tract the physical quantities from our phase field model. 
As discussed earlier, solutions to the extended gradient 
flow system (||) consist of grains in which V# = sepa- 
rated by a network grain boundaries in which V# 7^ 0. 
Although no rigorous proof of the inevitability of this 
solution structure exists, there are two plausibility argu- 
ments. First, once this structure develops, one can show 
that it will persist Jl2| . And second, such a solution can 
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be explicitly found for a flat stationary interface between 
two grains shown schematically in Fig. ^. This case cor- 
responds to the zeroth order of the formal asymptotic 
expansion as we show in what follows. 



hood. Also fllE 



c 



FIG. 2: A schematic representation of the solution for a flat 
static interface between two grains. £ is the coordinate normal 
to the interface. 



Fig. U is highly suggestive. We have a bicyrstalline 
system with a localized region where properties change 
(i.e. the grain boundary). Let us therefore consider the 
e — > limit of the gradient flow system (||). We ex- 
pect the width of the grain boundary region, defined by 
V# 7^ 0, to shrink to zero. We employ the standard 
method of matched asymptotics [ fi"5| which uses the small 
ratio of the interface width to its radius of curvature as 
an expansion parameter. In this limit, all functions and 
their derivatives vary slowly along the interface compared 
to their variation across the interface. The analysis can 
therefore be advanced by defining a coordinate normal 
to the interface and scaling it by e. The order parame- 
ters as functions of this new coordinate are termed the 
inner solution. They are expanded into a power series in 
e and the resulting hierarchy of equations is solved order 
by order, subject to a matching condition with the outer 
solution. This outer solution, obtained by setting e to 
zero, is valid far away from the interface. 

The mathematics of the sharp interface limit of our 
model is atypical because the fields r\ and 9 obey two 
different sets of equations Eqs. (g) hold in the grain 
boundary region which is defined as a strip S between 
two smooth non-intersecting curves r+ (t) and T e _ (t) (see 
Fig. |). Outside this strip V0 = and 



= a 2 V 2 ? 7 - /„ 



(3) 



Following [[l7| , let us adopt curvilinear orthogonal co- 
ordinate system {r(x, y, t), a(x,y,t)}. Let r(x,y,t) be 
the distance of the point (x,y) in S from ri. On T e _, 
coordinate a is the arc length. We introduce a scaled 
variable £ = r/e, and expand r\ and 8 in a formal power 
series in e 

r](C,(T,t) = T) ((,(T,i) +er] l (C,a, t) + ... , (4a) 
e(C,a,t) = 9 (C,a,t)+e9 1 (C,a,t) + ... , (4b) 

This expansion is valid in S and its immediate neighbor- 



V = 


„ 1 d 


. d 


(5a) 


V 2 = 


i d 2 

e 2 d( 2 ' 


^£ +*(!), 


(5b) 


d 
dt 


v d 


1-0(1), 


(5c) 



where n is the unit vector normal to T e _ pointing into S (C 
increases in the direction of h) and t _L n is the unit vector 
along the lines of constant r. The curvature k — V 2 r is 
positive when h points away from the center of curvature 
of rf_. Normal velocity v — —dr/dt is positive when the 
interface moves in the direction of n. This configuration 
is shown in Fig. |3|. 



ve = o 



r + e 



\r(x,y,t) 

n a 




ve = o 

FIG. 3: The region S between the curves Fl and T e + in which 
V(9 / 0. 

In order to proceed we must fix the scaling of a, s, 
t v and Tg in the sharp interface limit. We select the 
scaling for which 1) a flat interface does not move, and 
2) all terms in the free energy Eq. (jl|) scale with the same 
power of e. We also assume that the mobility functions 
P and Q are independent of e in the sharp interface limit. 
These conditions fix 



a 



e a. 



es, 



e 2 f„ 



T() 



e 2 fe. 



(6) 



The final ingredient of the formal asymptotic analysis 
is the matching of the inner and outer solutions. The 
inner solution is found separately in the grain boundary 
S (V0 7^ 0) and in its immediate neighborhood (V0 = 0). 
The inner solution valid in S is denoted by T](C,, a, t) and 
6{Q,u,t) without superscripts. The piece of the inner 
solution valid in the interior of a grain just outside of S 
is denoted by T]^{Q,a,t) and #^(C, <r,t). It is matched 
with the outer solution rj — 1 and — const for £ — * oo. 
We must also match the two pieces of the inner solution 
and their derivatives at the boundaries Tj. of the strip S 
which are at £ = £±. 



IV. ZEROTH ORDER SOLUTIONS: 
INTERFACE WIDTH AND ENERGY 

Now, having detailed the formal method of asymptotic 
expansion, we proceed to examine the results of this ex- 
pansion, term by term. We begin, naturally enough, with 
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the zeroth order. The ultimate result of the zeroth order 
calculation with be a determination of (i) the interface 
width (ii) the surface energy and (iii) the value of the 
order parameter 77 in grain boundary. 

Without loss of generality we can shift the origin of £ 
so that it is in the middle of S, C± = ±(Co+eCi + - • • )• We 
first look at the grain interior region C > £ + . Substituting 
the scaling ansatz Eq. (^) and the e-expansion (^) into 
(||) and using (||) we obtain 



= 5 2 (^ } )" 



9« 



)' = 0. 



(7) 



These solutions must be matched with the outer solution 
r/ = 1 and 8 = 8 + , as C — > +oo. Since the zeroth order 
functions 770 arL< i $o are independent of a, the coordinate 
along the interface, and time t, they describe a flat sta- 
tionary interface. They should therefore be symmetric 
with respect to the center of the boundary C = 0. We 
can thus set 8+ = Ad/2, half the total misorientation 
(which is in general a function of time) . Using the fact 
that /(l) = we arrive at 



A8 



(8) 



After similar manipulations we obtain the equations 
valid in S for Q € [-(+, C+] 



= -Um)- 2^o)(0o) 2 -^o)0o,(9a) 



= [/^oK + ^o)]', 



(9b) 



Note that we assumed that 8' > in 5. This assumption 
proves to be unrestrictive since all measurable quantities 
depend on the square of this derivative. 

To obtain the boundary conditions for 770 and 8q we 
employ two ideas. First, due to the aforementioned sym- 
metry, 770 is even and #0 is odd in C- Second, the re- 
quirement that 77, 8 and their derivatives be continuous 
at ( = C+ can be shown to lead to the continuity of all 
terms in the e-expansion (0) at C = Co- We thus obtain 



^(0) = 0, 77^(Co) 
9 (0) = , 8 (Co) = 



\/2/(77 max ) 



(10) 



A8_ 
~2~ ; 



9 (Co)=0, (11) 



where 77, nax = 770 (Co)- Let us also introduce 7y m i n = 770(0). 
This designation reflects our assumption that r/' > in 
C G [0, Co] so that 77 m i n < 77 < ?7 max . We can prove that, 
in particular cases, this assumption is indeed justified. 

Using the last boundary condition we integrate ( pb[ ) to 
obtain 



~ gfanax) - g(m) 

Hvo) 



(12) 



Upon substitution of this expression into Eq. ( |9a| ) we 
discover that r)' is an integrating factor. Using the second 



condition in Eq. (|T^) we obtain 



2f-s 



~2 (9(Vn 



9f 



-,1/2 



(13) 



The first condition in Eq. (|10|) furnishes the relation be- 
tween T7 m in and 



/ \ / \ , V 2 f(Vmin)h(T] mi n) , . 

liVmax) = giVrmn) H ; • (14) 



Armed with this condition we can obtain an equation for 
?7min via 



AO f Co f v ' 
— = / d(8' = as drj 



9(r/max) - 9 



vxmn h^/2f - s 2 (g(r) max ) - g) 2 /h 

(15) 



where we took advantage of the monotonicity of 770. Once 
Vmin is found we can calculate the width of the boundary 



Co = a 



drj 



'Vraln V 2 f ^ (S^max) ~ gf I h 

and the interfacial energy 



(16) 



7 = SA8 g(Vraax) + 25 / d V y/2f 

+ 2a [ d V sj2f - P(g(r) max ) - g^/h. (17) 



Remarkably, these formulas reduce the flat boundary 
problem to quadratures. We are able to compute 77 m i n 
and 77 max from Eqns. (|l4j) and (|l5|), the interface width 
from Eqn. (|l6|), and the surface energy from Eqn. (|l7|). 

As useful as these expressions are, insight can be better 
obtained from analytic expressions, since the above four 
equations must typically be solved numerically. Thus, let 
us examine the behavior of the zeroth order solution in 
certain limits in which an approximate solution can be 
found. The integral on the right hand side of Eq. ( |T5| ) 
can be calculated approximately when it is small. We 
identify two such situations. 



A. Small h approximation to the zeroth order 
solution 

To make contact with Ref. f| in which the |V#| 2 term 
in the free energy is set to zero, we consider a limit in 
which its coefficient vanishes h — > 0. Expanding Eq. ( fl4| ) 
in powers of h we obtain to lowest order 



? ?max - ry min « V J in = G « 1, (18) 

* ilr) 

where f mm = f(r] m i n ), etc. Thus, as h — > 0, the difference 
between 77 m i n and ?7 max vanishes, while they remain well 
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separated from 1. Therefore, /, g and h and their deriva- 
tives are regular at the undistinguished point r/ min < 1. 

Let us define an auxiliary parameter y via 77 = ?7 m i n + 
Gy and expand Eq. ([D]) in powers of G. We obtain 



A0_ 

2 



a g™ 



dy 



5# 



y/y(2-y) S 9n 



(19) 



One could in principle now invert (|T^) to solve for 7y m i n . 
The width of the boundary Co can be then found with 
the same accuracy 



7T d Vh min 
Co ~ — — — — ► 0. 



2 a 9. 



mm 



(20) 



The interfacial energy is, in this same approximation, 



7 » a 5 mm A<9 + 2d / dr/V2/ 



(21) 



For a special case examined in Ref. 



/ = !(i- 

77) 2 , g = r] 2 , and ft = 0, Eq. ( |l9[ ) is easily invertible. 
Indeed, the expressions for ?y m i n and 7 coincide with those 
of Ref. g. 

We finally remark that within this approximation, the 
behavior of 7 and Co m the limit of small misorientation 
AO depends on the properties of g and h near 77 = 1. To 
see why this is true, let us look at Eq. (|l^ ) in the limit of 
small AO. Unless g v is singular at some value of 77 other 
than 1 (unphysical) , small A6 implies small 1 — ?7 m i n since 
/(l) = 0. We shall see in the next subsection that it is 
true in general. This fact is not surprising since small 
angle boundaries can be thought of as arrays of distant 
dislocations. 



B. Small A9 approximation to the zeroth order 
solution 



All parameters being fixed, the right hand side of 
Eq. ( p"5| ) can be small if and only if 

P = f/max - »?mm < 1- (22) 

Since /(l) = 0, it is clear from Eq. that 

A = 1 - Vmin < 1. (23) 

Also, since we chose f(l) = we can approximate it by 
a parabola f(rj) « — 77) 2 near 77 = 1. Then 



P% 



X 



(24) 



If the behavior of g and h near 77 = 1 is known, we can 
obtain a complete solution to the zeroth order problem 
in this limit. Instead let us focus on the scaling of the 
interface width Co an< i the interfacial energy 7 in the 
A8 — > limit. This scaling can be deduced without a 
complete solution of the zeroth order problem. 



Suppose that near 77 = 1 (0 < A C 1) 
gf n - X 13 , h min - A 2 ". 



Then from (24) we obtain 

p ~ X 1+ "~P. 



(25) 



(26) 



Note that since p < X by definition, the scaling exponent 
in Eq. ( p^ ) must be greater than or equal to 1. This 
means that in this limit, the zeroth order solution exists 
only when u) > /3. 

To determine the scaling of the right hand side of 
Eq. (|l5| ) we define a finite integration variable 7/ S [0, 1] 
via 77 - 77 min = py. Then 



di] 

d-9 



P 



X 



2/- 



-2 (S(f7max) - <?) 2 



(27a) 
(27b) 

(27c) 



- [(2/A - 2f min h m " 



s'(g - g min ) 21 



+ 25 ( 5 -.g min )V2/ min /i min 



We used the fact that u > (3 to est ablis h that the 
second term in the square brackets of (27c) dominates 
(when cj = /3 all three terms in these square brackets are 
equally important). Substituting scaling relations J27| ) 
into Eq. @ we obtain 



(28) 



We are now ready to calculate the behavior of Co and 7 
in the AO — ► limit. Substituting the scaling relations 
(E7I) into (HI) we obtain 



Co ~ Ai^-« ~ (A^)TO. 



(29) 



The interfacial energy consists of three pieces 7 = 71 + 
72 + 73 which scale differently with AO. We list them 
separately 



71 = 2d 



^77 V / 2/-S 2 ( 5 (77 max )-ff) 2 /^(30a) 



A 2+§(u-/3) 



72 = 2d / dr?V2/ ~ A ~ 

73 = sA^^max) ~ 3 min A0. 



(30b) 
(30c) 



We can make several observations. First, since w > /3, 
72 always dominates 71 in the A0 — > limit. Second, 
when /3 < —1 we can integrate ( p5^ to obtain g mm ~ 
A 1+/3 . When /3 = —1, we obtain in a similar fashion 
g mm ~ In A. For f3 > — 1 the behavior of g near 77 = 1 
is arbitrary. Thus, the scaling of the interfacial energy 
(whether it is dominated by 72 or 73) in the limit of the 
small misorientation can be freely controlled by adjusting 
the behavior of g and h near 77 = 1. 
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V. FIRST ORDER SOLUTIONS: INTERFACE 
MOBILITY 

Having completed our analysis of the zeroth order in 
the asymptotic expansion, we now continue on to first 
order. This order in the expansion will yield the velocity 
of the interface as a function of geometry, mobilities, and 
surface energy. From classical as well as order parameter 
models, we expect this motion to be by curvature p8| , 
and, as we show below, this is indeed the case. 

To begin, in order for us to establish the matching con- 
ditions between the two pieces of the first order solution, 
let us look at the equation for rj{ valid in C > C+ 

-KQ + « 2 «](^)' = 5?{^ ] )" ~ Vi ] fvv(Vo } )- (31) 
Progress can be made by noticing that r)f* — (jy^')' is 



can take advantage of this fact by multiplying Eq. ( |3l| ) 
by T]q^ and integrating over the grain interior. Using the 
matching condition with the outer solution which states 
that all derivatives vanish as £ — > oo, we obtain 



a 2 



c+ 



c=c+ 

dCKQ + fi 2 K ][(^ l) )f . (32) 



a solution of (31) with the left hand side set to 0. We 



As we mentioned before, all orders of the e-expansion 
and their derivatives must be continuous at C+ or, equiv- 
alently, at (q. Consequently, we may drop the (i) su- 
perscript from the expression on the left hand side of 
Eq. j3^). This condition will suffice to set rji at ( = ( . 
We will not need the boundary condition for B\ at ( = (q. 

Turning our attention to the boundary region £ £ 
[0, Co] we write down the first order equations 



- vt v Q ?7q = c?{rl{ + kt)' ) - f vn m - ^mis'o) 2 - K Q'oO'i - sg m m^o ~ s 9v e 'n ( 33a ) 
VT S rjlP9' — sng + nh6' + [h6[ + h v rji9' + sg n rji]' , (33b) 



where we used the fact that V • [-r^fr] = n. The couplings 



to obtain a solvability condition which determines the 
velocity of the interface v. Multiplying Eq. ([36]) by rj' 
and integrating over [0, £o] we obtain 



/, <?, h, and their derivatives, and Q, P in Eqs. (33) are 
evaluated at the zeroth order solution. 

We can fix the boundary conditions at £ = for the 
first order functions by noticing that rjx is odd while Q x a 2 \rjtfj x - ^o^iIq = ^ ( 7 7i* > )' ( 7 7o* > )' _ V^iVo^)" 

is even in £. Therefore 



C=Co 



r ?1 (0) = 0, 0j(O)=O. 



(34) 



Using the second condition we integrate (33b) to obtain 



Hvo) Q'i = -[ h nM Oo + sg r ,(vo)} Vi 
rC 

d(' (vfg P{r}o,fi9' )r)l6' - sk g(i] ) - k h(r/ ) 6' ) . 



(35) 



Upon substitution of this expression into (33a) we obtain 
an equation for 771 of the form 



C[rn] = S?rf{ + Cfo, o )»7i = D( m , 6 ), (36) 



where 



D(r) Q ,6 ) = -r] (vf v Q + k) 



h 



d(' [vf e p nlo'o + «(K + sg)]. (37) 



The exact form of C is unimportant since we can show by 
direct substitution that (as in a conventional asymptotic 
expansion problem) C[rj' Q ] — 0. This fact can be utilized 



„(0\>i2 



Co 



Co 



Co 



d( [vf n Q + a 2 K](r]' ) 2 



dCe [vfePvX + <h6 + sg)]. (38) 



Solving for the interface velocity v and using Eqs. (|12|), 
@ and (0) yields 



-kjM 



(39) 



As expected, v is proportional to both the curvature and 
the energy of the interface. The mobility is given by 



1 

M 



dCQ(vo,nO'M) 2 + f e I d(P( Vo ,n6' )ri 2 (6 ) 

(40) 



We dropped the superscript (i) since it is clear that solu- 
tions which are valid in the grain interior should be used 
in @ for C > Co- 

Just as in a conventional formal asymptotic analysis, 
we obtained the normal velocity of the interface without 
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having to solve the first order equations. This is a general 
feature of analyses of this kind and allows one to express 
the mobility of the interface in terms of the properties of 
a stationary interface. 

To better understand the behavior of the interface mo- 
bility M 1 let us apply the approximations of sections 
f\TS and ffVBl to Eq. pi). 



A. Mobility in the small h limit 



structure. In the sharp interface limit, the solution con- 
sists of a set of regions Gi of spatially uniform 9 (grains) 
separated by narrow (or order e) grain boundaries. To 
calculate the rotation rate fl; of grain i, we integrate 
Eq. (2b) over the grain interior Gi. We obtain 



r<A / dAr) 2 P(r),0) 



s (t> da g{rj) h ■ 



V6> 



s$>iS(^ax)sign(A%). (44) 



To illustrate the importance of the |V6*| 2 term in the 
free energy for the motion of boundaries, let us again 
consider the limit in which its coefficient h vanishes. The 
second term in Eq. ( |4C| ) dominates in this limit. Assum- 
ing for the sake of the argument that P does not depend 
on V0, we obtain 



1 IT Tg a P 



'/mi 



(41) 



The mobility of the grain boundary therefore vanishes as 
h 1 / 2 in this limit in support of our claim that the \V9\ 2 
term is required for migration of boundaries. 



B. Mobility in the small AO limit 

It is instructive to trace the behavior of the interface 
mobility in the limit of the vanishing misorientation. For 
simplicity we assume here that the mobility functions P 
and Q are regular and assume a non-zero value at rj = 1 
and V0 = 0. Using @ we obtain the scaling of the 
three pieces which make up the right hand side of (Eoh 



Co 



I \2 



dCQWo) 



Co 



Co 



dCPvUO'of 



A 2+ — /5 , 
A 2 +^, 



(42a) 
(42b) 

(42c) 



Eq. ([28|) can now be used to determine the behavior of 
M in the A9 — > limit. We obtain 



1 

M 



UAO) 
\(A9) 



2 + ui-P 
t-t£±0 



P>0, 
0<O. 



(43) 



To illustrate, if g and h are regular at rj = 1 so that 
u) = f3 = our analysis yields 7 ~ A9 while M~ l ~ 
(A9) 2 so that the interface velocity diverges in the limit 
of vanishing misorientation v ~ (Aff)^ 1 . 



VI. GRAIN ROTATION 

As we discussed above, this model also allows for the 
grains to rotate. Let us consider a fully developed grain 



Here Qj is the length of the common boundary between 
grains i and j. The summation is over the neighboring 
grains j of orientation 9j . In integrating by parts we used 
the continuity of b = j^|y at the edge of the grain . In 
the grain boundary S, b is a unit vector in the direction 
of increasing 9. Thus, away from a triple junction, b = 
hij sign(A6>ij ) on both edges T± of the boundary between 
grains i and j. The unit normal n,j points from grain i 
into grain j. The periodicity of 9 must be taken into 
account to calculate sign(A6>y ). 

Recall that ?y is exponentially close to 1 in all of Gi 
except a narrow (of order e) strip near the boundary. 
Therefore, if r\ = 1 is neither a zero nor a singularity 
of P, the integral on the left hand side of (|44|) is ap- 
proximately equal to P(l, 0), where Ai is the area of 
the ?-th grain. This choice of P is problematic, however, 
since in the sharp interface limit ~ s/tq ~ 1/e. We 
therefore reach a remarkable conclusion that when P is 
regular at r\ = 1, grain rotation dominates grain bound- 
ary migration in the sharp interface limit. The reason for 
this unexpected result is that 89 jdt is continuous across 
the edge of the grain boundary. Inside the grain bound- 
ary, this time derivative must scale with the inverse of 
the grain boundary width 1/e. The time rate of change 
of the orientation 9 in the grain interior must therefore 
have the same scaling. 

Another way of interpreting the divergence of the rota- 
tion rate in the sharp interface limit is to consider what 
happens in the interior of the grain. Since the orienta- 
tion order parameter is constrained to remain uniform in 
space, it may be thought of as obeying a diffusion equa- 
tion with infinite diffusivity. It is therefore not surprising 
that for a generic choice of the mobility function P, the 
rotation rate is of the same order as the rate of change 
of 9 in the grain boundaries (fast). The resolution of this 
paradox is in fact satisfyingly physical. When material 
is nearly a perfect crystal (rj close to 1), i.e. there are few 
defects, one should expect the rate at which the order pa- 
rameters change to vanish. This can be accomplished by 
letting the mobility functions P and Q diverge at r\ = 1. 

In fact, that for a particular choice of weakly singular 
P(r],0) ~ — ln(l — rj) the grain rotation rate no longer 
diverges in the sharp interface limit. The integral on the 
left hand side of (44) may be calculated approximately 
in the sharp interface limit. We will need an expression 

(i) 

for ?7q far away from the boundary where it is close to 
1. Using the approximation for / sa |/ w (l)(l — if) 2 , we 
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obtain 



CO 



exp --^—i C 



(45) 



Therefore, focusing on the e-scaling in the sharp interface 
limit, we deduce that 



J dA ?7 2 ln(l -77) - LeJ (d( 



5 



(46) 



where L is some macroscopic length of the order of the 
grain size. Therefore Qi ~ se/rg ~ 1 no longer diverges 
in the sharp interface limit. Another important conse- 
quence of this argument is that since the perimeter of a 
grain is proportional to L, the rotation rate Q is inversely 
proportional to L 2 



(47) 



This prediction is independent of the choice of other the 
model functions. It is inconsistent with the heuristic 
derivations jn], [l9| of the rotation rate due to the dif- 
fusion of atoms along the grain boundary. These stud- 
ies obtain a 1/L 3 or 1/i 4 scaling of fi depending on the 
mechanism. In a separate study, Martin |9J assumed that 
rotation is caused by viscous motion of dislocations and 
obtained a rotation rate which was independent of L. 



VII. APPROXIMATE MODEL OF A SINGLE 
CIRCULAR GRAIN 

To illustrate the predictions of the sharp interface limit 
calculation of the preceding sections, we consider a cir- 
cular grain of radius R and orientation 9 embedded in an 
immovable matrix of orientation 0. To make analytical 
progress we choose Q = 1, P = — ln(l — 77), / = ^(1 — t?) 2 , 
g = —2rj — 2 ln(l — 77), h = 1 to ensure finite rotation rate 
in the sharp interface limit and Read-Shockley [^0| be- 
havior low angle boundaries. We also restrict ourselves 
to the small 9 approximation. In this limit we can carry 
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FIG. 4: Evolution of the squared radius R 2 (t) of the circular 
grain for 02 = 1 and three values of 01. 
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FIG. 5: Evolution of the orientation 8(t) of the circular grain 
for the same three values of 01. 



out the expansion in detail to obtain 



1 - In — 

a 



7r &6 1 
^ ~ 4 VT' M " 25 2 



r v st 



(48) 



Applying the motion by curvature result (|39| ) we obtain 

2RR&ailn6, (49) 

where ai = 4a 2 /f n . The expression for the rotation rate 
with the radius of the grain R and misorientation 9 can 
be obtained via the arguments of this section above. We 
obtain 



»2 



In 9 



(50) 



where a>2 — Qsa/fe- Solutions to these equations for 
R(0) = 1 and 0(0) = 0.1 are given in Figs. | and | The 
ratio 01/02 controls the behavior of the solution. When 
this ratio is small, rotation dominates the dynamics so 
that the radius of the grain is not significantly reduced 
by the time the grain rotates into alignment with the 
matrix. On the other hand, when this ratio is large, the 
evolution of the radius squared of the grain is almost 
linear in time as in the case of the motion by curvature. 



VIII. DISCUSSION 

In this paper we analyze a modified version of the phase 
field model of KWC |J . This model is constructed to de- 
scribe rotation of crystalline grains coupled to the motion 
of grain boundaries. The order parameter 9 reflects the 
local crystal orientation, whereas rj represents local crys- 
talline order. The Ginzburg-Landau free energy depends 
only on V# and is therefore invariant under rotations. In- 
clusion of the non-analytic \V6\ term into the free energy 
results in singular gradient flow equations. However, this 
singularity can be dealt with in a systematic way. 

Quite generally, solutions to the model represent a col- 
lection of regions of uniform 9 — grains — connected by 
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narrow (of order e) internal layers — grain boundaries. We 
are able to calculate the velocity of the boundaries in the 
limit of vanishing interface thickness, and find that it is 
proportional to the product of surface energy, curvature 
of the interface, and a mobility which depends on model 
parameters. The behavior of the intcrfacial width, en- 
ergy and mobility in the limit of the small misorientation 
is controlled by the behavior of the model couplings near 
7] = 1. 

We calculate the rate of grain rotation in the sharp in- 
terface limit and find that it diverges unless the mobility 
function P is singular at r\ = 1. For a logarithmic choice 
of this singularity, the rate of the grain rotation is finite 
and non-zero in the sharp interface limit. We explain 
this mathematical requirement by noting that the singu- 
lar term in the free energy results in infinitely fast dif- 
fusion in the interior of a grain. Therefore, 77-dependent 
mobility P must compensate for that fact. This loga- 
rithmically singular P leads to the conclusion that the 
rotation rate of a grain scales with the inverse of its area. 
This is a robust prediction of our model independent of 
the choice of all other couplings. 



For a plausible choice of model functions, motivated 
by the physics of low angle grain boundaries, we derive 
and solve equations describing a circular grain embedded 
in a matrix. We find that, as expected, when the scaled 
coefficient s of the |V0| term in the free energy is large, 
rotation is fast so that the radius of the grain does not 
change much by the time the grain rotates into coinci- 
dence with the matrix. While when s is small, rotation 
becomes important only when the radius is significantly 
reduced. 

We conclude by remarking that the model may be read- 
ily generalized in a variety of ways. We may extend the 
model to three dimensions, by constructing an appropri- 
ate tensor order parameter which reflects the symmetries 
of the lattice. Anisotropy may be included by allowing 
the coefficient of the |V?7| 2 term to depend on V0. Alter- 
natively, the mobility functions may be made anisotropic 
to yield kinetics which depend on the orientation of the 
boundary region. Overall, this model will provide a foun- 
dation for a physical, yet still relatively mathematically 
simple, model of grain boundary evolution and grain ro- 
tation. 
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